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Abstract 

We prove global in time existence of solutions of the Euler compressible equations 
for a Van der Waals gas when the density is small enough in H m , for m large enough. 
To do so, we introduce a specific symmetrisation allowing areas of null density. Next, 
we make estimates in H m , using for some terms the estimates done by M. Grassin, 
who proved the same theorem in the easier case of a perfect polytropic gas. We treat 
the remaining terms separately, due to their non-linearity. 
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1 Introduction 

We are interested in the Cauchy problem for Euler compressible equations, describing the 
evolution of a gas whose thermodynamical and kinetic properties are known at time t = 0. 

More specifically, we are concerned with the life span of smooth solutions. Various 
authors, in particular T. C. Sideris [20, 21], T. Makino, S. Ukai &; S. Kawashima [13], J.-Y. 
Chemin [2, 3] have given criteria for mathematical explosion. We know also that there 
exist global in time solutions for well chosen initial data. Li Ta Tsien [12], D. Serre [19], M. 
Grassin [7] prove for example the global in time existence of regular solutions under some 
hypotheses of "expansivity". 

Most of these results were obtained within the framework of Perfect Polytropic Gases. 
A natural question is to determine whether these results extend to more realistic gases, 
following for example the Van der Waals law. This law takes into account the volume of 
molecules, which is important in physical situations like explosions or implosions. In such 
limits, the gas is highly compressed and the Van der Waals law fits better with the behaviour 
of real gases than the Perfect polytropic gases in such limits. The Van der Waals law is 
also used to modelise dusty gases, seen as perfect gases with dust pollution [8, 16, 22, 24]. 
This law is given by the relationship 

p{v -b) = mT, 

where p is the pressure, v the massic volume, T the temperature and b, 9\ are given con- 
stants. 
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The addition of the covolume b, representing the compressibility limit of the molecules 
in the gas, modifies non-trivially the analysis of the Euler equations. However, we are going 
to show the global in time existence of regular solutions thus generalising a theorem of 
M. Grassin [7]. 

Theorem 1.1. Let m > 1 + d/2. Let (po,uq,sq) be the initial conditions for the Cauchy 
problem associated to the Euler compressible equations (2.1) for a Van der Waals gas with 
constant c v and c v > 0. Let us assume ^ po ^ 1/6. Then we can define 70 = 1 + 
furthermore there exists £q > such that if 
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(HI) ||(7ro,so)|| H -( K d) < £ o, where ir = (jOfe) 2 exp(^), 

(H2) the initial speed uq belongs to the space I = {z:R + xM d 4 R d ; Dz G L°°, D 2 z € 
H m_1 }, 

(H3) there exists 5 > such that for all x £ M. d , dist(Spec(Dito)(x), R_) ^ 5, 
(H4) the initial density po and the initial entropy so have compact support, 
then the problem 

{ d t u+{u-V)u = on M+xIE^, 

1 u(0,x) = u Q (x) onR d . ^ ' ' 

admits a global classical solution. If furthermore 70 = ^j- with v € N and v ^ 2, or if 70 
and m satisfy v = ^ m > 1 + \, then there exists a global classical solution (p,u,s) 

to the Euler compressible equations (2.1) satisfying 

(^-JL-^J 1 , u -u,s) G^°(lR + ;H m (M d ;]R d+2 ))n^ 1 (lR + ;H m " 1 (M d ;lR d+2 )). 

To do so, we have first to extend to Van der Waals gases a symmetrisation obtained 
by Makino, Ukai & Kawashima [13], which allows null density areas. Next we will derive 
energy estimates in H m . 

Since the Van der Waals gases have a behaviour close to perfect polytropic gases for 
weak densities, this result is not much surprising. However, the non linear terms introduced 
by the Van der Waals law have to be treated carefully. 

In Section 2, we describe the thermodynamic properties of a compressible gas and we 
state some important properties such as the Friedrichs' symmetrisation. In Section 3, we 
give the detailed proof of this result and in Section 4 we give the proofs of some technical 
lemmas used in Section 3. 



2 Thermodynamic and first properties 
2.1 Conservation law 

Compressible fluid dynamics, without viscosity or heat transfer, is described by the Euler 
equations, which are made of the conservation of mass, of momentum and of energy (see 
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[5, chap. 2]) : 

d t p + dxv(pu) = , 

d t q + div (pu <g> u) + Vp = , (2.1) 
d t E + div((E + p)u) = 0, 

where p is the mass of the fluid per unit of volume, q = pu is the momentum per unit 
of volume and E = ^pu 2 + pe is the total energy per unit of volume, sum of the kinetic 
energy and internal energy. This is a system of (d + 2) equations and (d + 3) unknowns: 
the density p G M+, the speed u 6 R d , the internal energy e G M and the pressure pGl. 
In order to complete this system, we have to add a state law, for example an incomplete 
state law, also called pressure law (p, e) i-> p(p, e). 

Definition 2.1. We call Euler compressible equations the system made of the conservation 
laws (2.1) and of an incomplete state law p = p(p, e). 

A simplified model is often considered, conserving only the conservation of mass and 
momentum, assuming that the fluid is isentropic. This simplified system is 

| d t p + div(pu) = , (22 . 
I d t q + div(/9n <g> u) + Vp = , { ' ' 

and the state law is a given function p i-> p(p). 

Definition 2.2. We call isentropic Euler equations the system (2.2) with a given state law 
p = p(p) such that 

de 

^ dv 

where v = 1/p is the specific volume, T is the temperature, s is the specific entropy and 
/ = e — Ts is the specific free energy, assuming we are given a complete equation of state 
(v, s) ' y e(v, s). 

The thermodynamical variables v,s,e,T,p must satisfy some relations described in Sec- 
tion 2.2. The Euler equations are a system of first order conservation law, whose study is 
developed in particular in the books [4, 17, 18]. 



2.2 State law 

The state law has a strong influence on the mathematical analysis of the compressible Euler 
equations. The state law of "real" gases can reveal particular behaviour and introduce 
existence and/or uniqueness troubles which do not appear for perfect gases, see [14]. We 
describe below the physical principles a state law has to satisfy. 



2. 2. A Definitions 

We consider a fluid, whose internal energy is a regular function of its specific volume 1 
v = 1/p and of its specific entropy s. This means that the gas is entitled with a complete 
state law, or energy law e = e(v,s). The fundamental thermodynamic principle is then 

de = -pdv + Tds (2.3) 

1 specific is a synonym of massic 
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where p is the pressure and T the temperature of the gas. Consequently, the pressure p 
and the temperature T can be defined as 



P 



de 
dv 



de 
ds 



(2.4) 



where the notation | precises the variable maintained constant in the partial derivation. 

The greater order derivatives of e have also an important role; we introduce the following 
adimensional quantities: 



v dp 
p dv 



r = -- — 

T dv 



pv dT 
T 2 ~ds~ 





d 3 e 




V 


dv 3 


s 


2 


d 2 e 






dv 2 


s 



(2.5) 



The coefficient 7 is called the adiabatic exponent, and T is the Griineisen coefficient. The 
quantities 7, 5, T and ^ characterise the geometrical properties of the isentropic curves in 
the (v,p) plan (see [14]). They can be expressed in function of e through the relationships: 



7 



v d 2 e 
p dv 2 



v d 2 e 
Tdsdv ' 



pv d 2 e 



We also introduce the calorific capacity at constant volume c v and the calorific capacity 
at constant pressure c p by 



de 
dT 



d 2 e 




W 1 


V 



ds 

Cp= dT 



(2.6) 



These two quantities are linked with ^ and with 7, 5, F through 

pv pv 7 



5c v 



T 7 <5 - r 2 

7<5 



(2.7) 



The quantity 7* = can besides be expressed as 7* = g_ r $ ■ It is not equal to 7 in 
the general case, but for an ideal gas we have 5 = V = 7 — 1, so that 7* = 7. 



2.2.B Thermodynamical constraints 

It is very natural to assume v ^ 0. We assume furthermore that the pressure p and 
the temperature T are positive, which imposes that e is a function increasing in T and 
decreasing in v. 

A classical thermodynamical hypothesis requires furthermore e to be a convex function 
of s and v, which means: 



j5 - r 2 ^ , 



7^0. 



In particular, 7^0 means that p increases with the density p = 1/v, which allows us to 
define the adiabatic sound speed by 



dp 

d~ P 



'7- • 
P 



(2.8) 
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Then, we show that ^ can be expressed in function of p and c through the expression 
y _ 1 9 (P C ) 

c dp s 

Furthermore, we require usually T > and > 0. The condition T > is not 
thermodynamically required but is satisfied for many gases and ensures that the isentropes 
do not cross each other in the (v,p) plan. The condition <£ > means that the isentropes 
are strictly convex in the (v,p) plan. 



2.2.C Van der Waals Gas 

Definition 2.3. A gas is said to follow the Van der Waals law, if it satisfies the following 
pressure law: 

p(v-b)=WT, (2.9) 

where v is the massic 2 volume and b is the covolume, representing the compressibility limit 
of the fluid, due to the volume of the molecules. 

The Van der Waals law is a modification of the perfect gas law, in which 6 = 0. In 
opposition to the perfect gas law, it takes into account the proper size of the molecules, 
which is important in some situations when the gas is strongly compressed. In this model, 
the density must be bounded and the maximal density is p ma x = \- 

The fundamental relationship (2.3) gives us the PDE: d v e + -^rd s e = 0. Thus, we 
introduce new variables w = (v — 6) _iH , a = (v — 6) _K exp(s) and e(w,a) = e(v,s). We 
obtain d w e = 0, so that e = £((v — 6) _9 *exp(s)) for any regular function £. 

If we assume furthermore that c v is constant, thanks to the definition of c v and (2.3), 



we get that 



d 2 e 



Cv ds 



hence aS" = - !)£' and £{a) = Ca x l c - which leads to: 



e = (v — b)~~ exp( — ) , p 



' c v c v v — b 

Some computations allow us finally to obtain 

v , . v ^ 70 + 1 t) 

7 = 7o t, T = 5 = (70 - 1) 



v — o v — b 2 v — b 

where 

7o = — + 1- (2-10) 

c v 

The conditions of Section 2.2.B are then satisfied for 70 > 1. 

Remark 2.4. 1. A perfect gas can be seen as a Van der Waals gas with b = 0. A perfect 
gas for which c v is constant is called polytropic. 

2. The Van der Waals law coincides with the dusty gas law [8, 9, 16, 22, 24]. In this 
model, we consider that the gas is perfect but polluted by dust particles that are 
equidistributed and have a non-negligible volume. 

Very often in the literature, the perfect polytropic gases are considered as a canonical 
example. However, their adequation with physical observations is not as good as for Van der 



2 also called specific 
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Waals gases, for example in explosion phenomena, or in the sonoluminescence phenomenon 
[1, 6, 11]. 

In the following, we only consider Van der Waals fluids with constant and strictly 
positive calorific capacity c v : 

c v >0, (2.11) 

which implies 70 := ^ + 1 > 1. 
2.3 Symmetrisation 

An important property of the Euler equation is their symmetrisability. 



2. 3. A General case, without vacuum 

dp 
dp 



If p > and & 



> then the system (2.1) can be written in the variables (p, p, s). Then, 

the system is almost symmetric, since it can be written matricially dfV+J2k Ak(V)dkV = 0, 
with V = (p, u T , s) T and 

( u-i pc 2 i T \ 

k y u-t j 

This matrix is almost symmetric since we obtain a symmetric matrix by multiplying it 
on the left by D := Diag p,...,p,lj. Consequently, we have the following 

Proposition 2.5. The system (2.1) is Friedrichs symmetrisable when (p,u,s) takes values 
in a compact subset ofW + x R d x R. 

Indeed, for such values of (p,u,s), DA(£,V) is symmetric with D symmetric definite 
positive. 

2.3.B Van der Waals gas 

We are not completely satisfied with the previous formulation as it does not authorise p 
to vanish, or even to tend to at infinity. For example, p cannot be taken into H m for 
m ^ 0. Makino et al. have introduced in [13] a symmetrisation for perfect polytropic gases 
allowing the null density areas. We generalise here their method to the case of Van der 
Waals gases. 

First, let us remind that for Van der Waals gases, we have 7 = and p = (70 — 

1) (i^p) exp (jh^- We now introduce the new variable 



7T 



70 / P \^o" 



7o-l \7o - 1, 

and we re-write the system (2.1) in the variables (ir, u, s). In order to do that, we first write 
the system (2.1) in (p, u, s) variables: 

dtp + u ■ Vp + pc 2 div u = , 

d t u + (u-V)u + ±Vp = 0, (2.12) 
dts + u ■ Vs = . 
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Since n = f(p) = Cp a , it is sufficient to multiply the first line by f'(p) = Cap a 1 to 
obtain an equation in ir: 

d t TT + u-Vir + Cap a - l pc 2 V -u = 0, 
5 t « + («.V)«+^ 5=T exp(^)V7r = 0, 

dts + u- Vs = . 

Besides, we know that c 2 = jz^j, C = ^J^i(^t) 270 and a = ~2^T- Jt remams to 
evaluate the coefficients: 

Cap pc = — 7r 

1 — bp 

70-1 7T 



2 1 - bp ' 



exp(s/(7 c 1 ,)) (70 - l) 1/7 ° i-"-i/7 

7T Q 



pCap 01 ' 1 " 1 - bp aC(T°-i)/(7oa) 

S \ 70 — 1 7T 

exp 1 



Thus, 



7 c„ / 2 1 - bp 



d t u + (u-V)u+ ^ T ^exp(^)V^ = 0, (2.13) 

c\s + u-Vs = . 

1 2 



Moreover, 1/(1 - bp) = 1 + b (-4^) 70 1 exp 7Tto-i . Therefore, denoting 



1 




,= ^>1, and b = b(^r) 70 -\ (2.14) 
7o-l V 4 7o y 

we obtain jjg^ = ^(l + be - ^ 70 ^:/!"^ 1 ), which is well defined for all ir ^ and in particular 
for p = 0, since we have assumed in (2.11) 70 > 1. The matrix associated to the system 
(2.13) writes now 

A(e, w )= 7o_i_^_ exp( _A_ K u . e7d 
\ 

It is once again "almost symmetric" in the sense that SA is symmetric, S being the diagonal 
definite positive matrix 

S = Diag (1, exp(-s/ {j c v )), . . . , exp(-s/ (70^)), l) . 

Furthermore, this symmetriser is independent from p and in particular is well defined and 
definite positive even when p or p vanishes. Finally, we have the proposition: 

Proposition 2.6. For a Van der Waals gas with constant c v and c v > 0, the system of Euler 
equations can be written for regular solutions as (2.13) which is Friedrichs symmetrisable 
for (p, u, s) € K, where K, is a compact subset o/R+ x R d x R. 
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Note that in the variables (n,u, s), the system (2.13) can be written 



d t it + u • Vtt + 22^1 ( 1 + b exp( 



70 c v ' 



d t u + {u.V)u+ ^ exp( _^ )(1 + 5 exp( _^ K -i )7rV 7r = 0, 

dts + u • Vs = . 



where and 6 are defined as in (2.14). 



(2.15) 



2.4 Local existence 

The symmetrisation of Proposition 2.6 is very useful to show local existence of regular 
solutions with vanishing density. 



Theorem 2.7. We consider a Van der Waals gas with constant c v such that 70 €]1,3]. Let 

9 po V 70 ~ 1)/2 ,70-1 

7r = 2 A / — exp(- s ), 

where p G tf l (R 3 ; [0, l/b[). We introduce also v = ^±±, b = b (^f}~ ■ We assume 

that (tto,uq, so) G H m (M. d ) for m > 1 + |. Then there exists T > anc? a unique solution 
(n,u,s) G ^([OjT] x IR 3 ) to i/je Cauchy problem for 

d t TT + u-VTT+ 3Scl(l + 6vr^ 1 )7rV-n = 0, 
5 t n + (n-V) U + 2^i e xp(^_)(l + 6exp(^)7r"- 1 )7rV7r = 0, (2.16) 

C?tS + m- Vs = , 



wif/i initial condition (ttq, uq, sq). Furthermore, 

(vr, «, s) € tf([0, T] ; H m (R d ; M d+2 ) ) n 'T 1 ( [0, T] ; H m_1 ( 



! ))- 



2.5 Positivity of the density 

For regular solutions, the positivity of the density is given by the following 

Proposition 2.8. Let (p,u,s) G ^([O, TexfxlR^) be a regular solution of the Cauchy 



problem (2.1) associated to the regular initial conditions (pq,uq,sq) G ^(W 1 ). If Vii G 
L°°([0,T] x R d ) for all T < T cx , p G L°°(M d ,M) and if p (x) ^ for all x G M. d , then for 
allte [0,T CX [, 

p(t,x) >0. 

Proof. We use the characteristics method to obtain an expression of the solution of the 
Cauchy problem (1.1). Assuming Vu is bounded, we obtain: 

p(t, x) = po(X(0; t, x)) exp ^— J div n(r, X(t; t, x)) dr^j >0, 

where X is a solution of the Cauchy problem 



which is global in time since Vit is bounded. 



X(t ;t ,x ) = x , 



□ 
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We also prove that for a Van der Waals gas with constant c v > 0, the variable tt = 

270 introduced in section 2.3.B remains non-negative if it is non- negative 

at initial time. This property implies in particular that if p$ < 1/6, then, as long as the 
regular solution exists, this property will be satisfied. 

Proposition 2.9. We consider a Van der Waals gas with constant c v . We denote tt = 

^i(^t) ^ ^=^T an ^ = H^) ^ • Let(ff,u,s)€^H[0,T cx [>cR d ) be a 
regular solution of (2.15) satisfying the initial conditions 

7r(0, x) = TTo(x) , u(0, x) = u (x) , s(0, x) = s (x) , 

with (tt ,uo,s ) G ^(R 4 *) and s G L°°. 7/divu G L°°([0,T] x R d ) for all T < T cx , 
vr G L°°(R d ,R) and < p (x) < 1/6 for all x G R d , then for all t G [0,T CX [ ; vr > 0. T/ten, 
we can define p, and we have: 

p(t,x) < 1/6. 

Proof. Let To < T cx . We introduce the Cauchy problem 

dtw + div (wu) = g(t, x, w) , (2-17) 
u>(0, x) = tw (x) , 



where 



s(t,x,ti;) = ( 1 - ^— ^(l + 6exp( g(t ' x) )^~ 1 ) ) «,div«(t,x) 



We can apply the Kruzkov theorem [10, 17]. Indeed, the hypotheses ensure that g(t, x, w) — 
wdiv (u(t, x)) = — ^^(1 + be~ s ^ 10Cv ^w v ~ 1 )w div u is uniformly bounded with respect to 
x G R rf when w is considered as a variable taking values in a compact set. Furthermore, 7To G 
L°°, so the regular solution 7f coincides with the entropy solution w\ of (2.17) associated to 
WQ,1 = ^o- 

Besides, the entropy solution u>2 of (2.17) associated to the initial condition u>o,2 = is 
the function constantly equal to 0. 

After Kruzkov Theorem too i ^ w o,2 implies twi ^ u>2 for all (i, x) G [0, To] x R rf , that is 
to say fr(t, s) > for all (t, x) G [0, T ] x R d . The formula: 

\ 




2 



1 + 6ex p(^;) 7r7D " 1 / 

allows us to conclude. □ 



3 Proof of Theorem 1.1 

In order to prove this theorem, we adapt M. Grassin's idea [7]. First, we look to the 
isentropic case, which allows us to simplify the estimates. For a Van der Waals gas, non- 
linear terms now appear in the estimate which we need to treat separately. 
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3.1 Isentropic case 

Let us consider first the isentropic case 



d t K + u-V-K + ^(1 + \m v - 1 )nV-u = 0, 
d t u+{u- V)n + 22 2 ll (l + &7r iy - 1 )^V7r = 0, 



(3.1) 

_. . , _ ,, 



with initial conditions 

7r(0, x) = vro(x) , u(0,x) = uo(x) . (3-2) 

which is technically simpler than the general case, but provides estimates very useful in 
order to treat the general case. 

We consider besides the problem 

f d t u + {u-V)u = on K + xR d , 

1 u(0,x) = uq(x) onR d , [ ' 

obtained by neglecting it in (3.1). After [7, Lemme 3.1 and Prop. 3.1] we have the following 
preliminary result 

Proposition 3.1. Under hypotheses (H2) and (H3), the problem (3.3) admits a global 
regular solution u satisfying 

1. Bu (t, x) = ^ + f^0, for all x G R d and t G R +; 



(i+t) 

2. ' 



3. 



B l u(t,-) s$ Ki{\ +t) d / 2 -('+ 1 ) forle [2,M + 1], 

L 2 

D 2 u(t,-)|| Loo < C(l + t)~ 3 , 



where K : R + x R d -> .M^W, II-^IIl^r+xI^) ^ C anc? ^ / or ' G rrx -h 1JJ &em# 
non-negative constants not depending on m,d,S, \\uo\\ x . 

3.1. A Local uniqueness 

Proposition 3.2. Let Uq = (ttq,uo) t G H m (R d ) and Uq be two initial data for (3.1). Let 
U = (ir,u) T , U be two corresponding solutions, defined for ^ t ^ To- We assume that 

< oo. Let xq G M d and R ^ 0. We denote 



DU 



L°°([0,To]xl 



M = SU p{(^-i|7r|(l + 67r^ 1 ) + | U |)(t,x),(t,x) G [0, T ] x B(x , R)} , (3.4) 
C T = {(t,x) G [0,7] x B(x ,R-Mt)} for T G [0, 7i] , (3.5) 

w/iere 7\ = min(T , r?/M). //17 = U on B(x ,R) then U = U onC Tl . 

The proof of this proposition is classical for hyperbolic systems (see for example [15]), 
the constant M being the maximal propagation speed. 

Proof. Let Uq = (ttq, uq), Uq = (ttq, uq) be two initial data for (3.1) such that Uq G H m . Let 
U, U be the solutions of the associated Cauchy problems. We assume that these solutions 
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are defined on [0, T ] with T > 0. Let also x € M. d , 7] € R* + and M, C T be as in (3.4) and 
(3.5). Then we have: 

d t U + ^a 3 {U)d 3 U = 0, 

3 

where 



a 3 (U) 



7Q- 1 





(tt + bTr") 



V 
















Consequently, 3 t (f7 - 17) + £\ a 3 (U)d 3 (U -U) + {a 3 (U) - a 3 (U))d 3 U = 0. Then, we make 
the scalar product with ({7 — £/) and we integrate on Ct for T G [0, To]. We get 

2 



2.L 



St C7 - C7 



+EM (f/ ~ c7) ' a?(f/)(c/ ~ c7) ) 



-^2(U -U) ■ dj( aj (U))(U -U)dxdt = - I (U - U) ■ (a 3 (U) - a 3 (U))d 3 U dxdt . 

3 J ° T 

Then using the Stokes formula and noting that 8Ct = ({0} x B(xq, rj)) U ({T} x B{xq, tj — 
MT)) UC, we obtain 



B{x ,R-MT) 



u-u 



'B(xo,n) 



U-U 



(0, x) dx 



! 1 f 

(T,x)dx-- \ 

1 Jb\ 

+ ; 1 / C - u\ 2 + Y(U - U) ■ a 3 (U)(U - U)^-du 

= / J2(U - U) ■ d 3 (a 3 (U))(U -U)-(U -U)- (a 3 (U) - a 3 (U))d 3 U dxdt . 

JCt j 

Besides ||a,-aj(u)|| LOO < C||DC/|| Loo(Ct) (1 + \\M\l^(c T ))- Hence ' 

I (U-U)- d 3 {a 3 {U)){U -U)-(U-U)- (a 3 (U) - a 3 {U))d 3 U dxdt 

J CV 



^C||D[/|| l ^ (Ct) (1 + ||<-^ t) 
Furthermore, the choice of M implies 

_ 2 



JB(x ,R-Mt) 



L 



u-u 



U-U 



dxdt . 



+ £([/ - U) ■ a 3 (U)U -U)j^-da>0, 

j 



so finally we get the estimate 



2^ 



B(x ,R-MT) 



U-U 



(T, x) dx - - 

z Jb(x , v ) 



^C||Di7|| LOO(CT) (l + ||7r||^- 1 (CT 




U-U 
U-U 



J B(x ,R-Mt) 



(0, x) dx 

2 



dx dt 
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We conclude thanks to Gronwall lemma that 

U-U 



R-MT) 



2 \ r 2 

(T, x) dx ^ -e c ' T / Uq-Uq (x) dx 

2 JB(x ,R) 



where C = C||D17|| Loo(Ct) (1 + \\M\l~\c T )) ■ ° 



3.1.B Local Existence 

We construct a local solution of (3.1)-(3.2) such that the difference between this solution 
and (0,n) be in ^°{R + ;U m (R d ;R d+1 )) n ^ 1 {R + ; H" 1 " 1 ^; R d+1 )). The first step is the 
symmetrisation of the system, given by Proposition 2.6. This result allows us to use a 
general theorem (see Theorem 2.7) giving the local existence of solution. Let us define as 
above, 

TQ-l 7Q-1 



2 



where we assume ^ p < 1/6. Then we use the same method as M. Grassin [7] to prove 
that the system (3.1) admits a local in time solution, with initial condition no in the space 
X and not in a Sobolev space (in particular, no does not tend to at infinity). We use here 
the compactness of the support of pQ (hypothesis (H4)) and the finite propagation speed of 
the solutions for an hyperbolic system. More precisely, we assume that Supp(po) C B(0,R) 
for R > 0. Let 77 > and ip G ^ c °°(R d ; R + ) be such that ip = 1 on B(0, R + 2r/). We obtain 
a local in time solution (tt^, u v ) of the problem (3.1) with initial conditions (710, <puo) G H m 
for i G [0,T[. The Propositions 2.8 and 2.9 ensure that the condition ^ p < 1/6 is 
satisfied. 

Let e G]0, T[ , we introduce the maximal propagation speed M = supj 70 ^ 1 (|7r y | + |7r^|' y ) + 
K| ; t G [0, T - e] , x G We also introduce e' g]0, ^[ and T x = min(T - e, ^ - e') 
the time for which this construction is available. We finally obtain a solution (71", n) of 
(3.1)-(3.2) by denoting 

(7T^,n^) in K, 
(0, n) out of K, , 



(7T,n) 



where K. is the cone /C = {(t, x) ; ^ i ^ Ti , x G i?(0, ii + r/ + Mt)}. Then it is sufficient to 
show that the solutions can be glued smoothly along dlC. We use here the property of local 
uniqueness given below by Proposition 3.2. Let indeed xq G S(0, R + rj) be the sphere of 
radius R + rj of centre and E XQ = {(t, x) ; t G [0, Ti], x G B(xo, 7] — Mt)}. The choice of Ti 
implies in particular dK, C U^gsmR+^-Exo ( see Fig- !)• The initial conditions of (tt^,^) 
and (0,n) coincide on E XQ since the support of ttq is included in B(0,R). Proposition 3.2 
allows us then to claim that (ir^,^) = (0,n) on E XQ and consequently on a neighbourhood 
of dlC. 

3.1.C Energy Estimates 

When we compare to the Perfect Polytropic Gas case [7], we observe that the system (3.1) 
has an additional term, which will modify the estimates. 
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Figure 1: Gluing of the solutions along d)C. 



Let us denote U = (it,u — u), w = u — u and U = (0, u). We also introduce 
( u-i — u-i 



Aj(U) 





7o -1 



7T 



Uj — Uj 




^vr \ 




B(DU,U) 



^t-tt div u 
(w • V)u 



Uj - Uj j 



(3.7) 



so that the system (3.1) can be written 



d t U + Aj{U)djU = -B(BU , U) - ^UjdjU - F(BU , DU , U) (3i 

3=1 3=1 



where 



F(BU ,BU,U) = t diY ( ^ + U) 



with v = 



7o + 1 
7o-l 



(3.9) 



Observing the properties of u described in Proposition 3.1, we expect the terms D k U 

for k G [0, m\ to decrease with respect to time with a rate depending on k. Consequent 
we introduce 

m 



L 2 

y. 



Y k (t)= D k U(t,- 



L 2 



Z(t) = Y,(l + t) 9 "Yk(t), (3.10) 



fc=0 



with gk = k + c, in which c has to be chosen so that all the terms of Z have the same 
decreasing in time. In order to estimate Z, we apply the operator D fc to (3.8), we make the 
scalar product with D fc U and we integrate on M. d . The system (3.8) is different from the 
one considered by M. Grassin through the term F(DU ,T>U ,U) defined in (3.9). We use 
now [7, Prop. 3.2 and 3.3] to estimate the terms in common. We remind these results here: 
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Proposition 3.3. Let a G N d be a d-uplet of size k ^ (that is to say \a\ l = a\ + . . .+dd 
k). Let us denote, for U = (ir,u — u), 



R k 



J B k U ■ B k A i (U)djU j dx , 

J B k U ■ B k ^B(BU ,U) + Y^ UjdjU^j dx 



then there exists a constant C > depending only on m, d such that 



\R k \ ^ C(l + t)-^- d ' 2 Y£z , S k ^ C(l + t)^- 2 Y k Z - ^-Y h 



l+t 



k 1 



where 



and 



r = min(l-^,(y-l)(i). 
We have now to estimate - J Rd B k U • B k (F(BU , BU , U)). Let us denote b 



(3.11) 



2 u 



so that - J Rd D fc J7 • D fc (F(Df7 , BU , £/)) = L + J. A priori, J is easier to estimate than I. 
However the estimate of / is possible since the matrix ir u ( ^ ^ J is symmetric. 



a) Estimate of /. We show in Section 4.1 the following 
Lemma 3.4. With the notations introduced in Section 3.1 

2 

L 2 



|/| «S CIIttII^IIDC/Hlc 



B k U 



(3.12) 



b) Estimate of J. Here, we divide J in two parts: a first part J\ which contains only 
first order derivatives of u, and a second part Ji in which all the derivatives of u are at 
least of order 2. More precisely, 



Ji = -b[ V (d a 7r)(d a TT u )dwu, 



a6N d ||a| 1 =fe 
J 2 = J~Jl. 

For Ji, we use the first point of Proposition 3.1 giving the decreasing in time of Bu, and 
Lemma 4.7 giving the estimate: 



|d a (^)|| L2 < llvrll^ 1 



B k TT 



L 2 



for all eZ-uplet a G N d of size A;, that is to say satisfying \a\ 1 = Yl a i = k. 
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We obtain, using the Cauchy-Schwarz inequality: 

|Ji| ^6||a a 7r!| L2 ||a a K)|| L2 ||divH|| LO 



c 



l + t 



B k U 



L 2 



7T 



(3.13) 



For J2, we prove in Section 4.1 the following 



Lemma 3.5. With the notations introduced in Section 3.1, there exists a constant C > 
such that 



\J 2 \ ^ C(l + t) 
w/iere 4 = (-g 1 - f + - 1) - g fc - 2. 



L 2 



(3.14) 



c) Re-assembling of the estimates. Assembling the results of Proposition 3.3, the 
estimates (3.12)-(3.13)-(3.14) of / and J, and finally using Lemma 4.5, we obtain 

1 <jy* + !L±1yZ < C(l + i)^ + C"(l + t)-^~ 2 Y k Z + C(l + t )P+V+iK»-i) Y *Z'' 

2 cl£ 1 "I - t 

+ C(l + t)(^+i)("-i)-iy2^«'-i + C (l + i)(/3+i)(^-i)-9fc-2y fc ^ ; (3, 15 ) 
where we have denoted 

Here, we choose the constant c introduced in ^ = k + c in order to have (3 = and 
consequently a good decreasing in time. This means to require g\ + d/2 = and c = 
— 1 — d/2. Consequently, we have 

1, d 1 
9k = k---l. 



We introduce now 



a = l + <i/2 + r>l 



so that k + r = g^ + a. We can now divide by in (3.15), multiply by (1 + t) 9fc and 
summate on k to obtain a differential inequality in Z (defined in (3.10)): 



dZ 
~dt 



+ ^ z<c (W_Z_ + (1+t) ,-. z , +1 ) 



Then, we introduce ((t) = (1 + i) a exp (^ij and we deduce from the inequality just 
above 



dC 



C 



dt (l + t) 

Besides, ( 2 + ( u+1 < 2£(1 + ( u ) for v ^ 2. Therefore 



1 dC < C 



that is to say 



C(i + C) d * " (i + *) a ' 
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with f(x) = i In . By integration, we obtain /(£(*)) + ^(1 + 1)"^ 1 ) < /(C(0)) + 

As / is strictly increasing and one-by-one from to M*_, if /(C(0)) + C/(a — 1) 
belongs to the set on which f^ 1 is well-defined, we obtain 

C(t) < Z" 1 (/(C(0)) + C/(a - 1)) . 

But /(C(0)) + C/(a— 1) ^ is only possible if £(0) is small enough, since f(x) tends to — oo 
when x tends to 0, and ("(0) = exp (C) Z(0) = exp(C)||7ro|| Hm - The smallness condition is 

satisfied thanks to the hypothesis (HI) with < Eq < / _1 (j^zfj- 
d) Conclusion. We have obtained the following inequalities 

n(t) < (i + t)-»z<() 

<(T^ exp (rS) rl ( /(eCz(0)) + ^)- 

The L 2 norms of the derivatives of the local solution U consequently do not explode in 
finite time since t i— >■ jj^e exp do not explode in finite time. Let us assume that 

the regular solution exists to time T. Our estimates give us, for all t <G [0, T[, for Ct not 
depending on T, 

||(7r,u-u)(i)|| Hm ^C T . 

Since Supp7r C B(0,R), our construction is possible when the norm of (ir,u) is bounded 
in H m (M d ) and 

II (*.«)(*) II H«.(B(0,fl)) < ° T+ \H^\\nm(B(0,R)) < R T- 

We can associate to the constant Kt a time of existence T*(Kt) for the local in time 
solution. Let t\ G]0, T[ be such that t ^ T — T*(Kt)- Introducing the solution with initial 
condition (w(ti),u(ti)), we succeed in prolongating the solution up to time T, which finishes 
the proof. 

3.2 General Case 

3. 2. A Local in time Existence 

As in the isentropic case, we first seek to symmetrise the system. Let us denote 



7o ~ 1 ( P \ 270 = / 7o ~ 1 f P \ 2 exp ( 7o - 1 s 



7o V7o-l/ V 7o V 1_ W V 2 7o c v 

The system (2.1) can be written in variables (tt,u,s): 

7Q + 1 

e s/( l0 c v ) dt7T + e s/( 10 c v ) u . v?r + 20_l e s/( l0 c v ) 7Tdivu = _10z±l JrJO -i diyu ; 

I ~ 7Q + 1 



d t u + (u ■ V)u + ^^V/^^ttVtt = -^for™- 1 Vtt, ( 3 ' 16 ) 

(i + t)- e (a t s + u- vs) = o, 
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where b = b y~%^~J ~'° 1 • We introduce furthermore a parameter 6 to be determined so that 
(1 + s has a decreasing in time similar to the estimates obtained in the isentropic case. 

In order to obtain local existence of a solution, we construct a solution by following the 
same strategy as in the isentropical case and using once again a property of local uniqueness, 
given by Proposition 3.6. 



3.2.B Local in time Uniqueness 

We show here a similar result to the one obtained in Section 3.1. A in the isentropical case. 



Proposition 3.6. Let Uq = {ttq, uq, sq) t G H m (M. d ) and Uq be two initial data for (3.1). 
Let U = (it, u, s) t , U be the two corresponding solutions defined for ^ t ^ Tq. We assume 



that 



DU 



L°°([0,T ]xl 



< oo. Let xq G M. d and R ^ 0. We denote 



M = sup{e 2 ^(2^i| vr |(i + fe|^|- 1 ) + | n |)( t)a; ) ; (t )a; ) G [o,T ] x B(x ,R)}, (3.17) 



C T = {(t, x) € [0, T] x B(x , R - Mt)} for T G [0, Ti] , 

where T\ = min(To, rj/M). 

If Uq = Uq on B(xq, R) then U = U on Ct x ■ 



(3.18) 



Proof. Let Uq = (ttq, uq, so), Uq = (ttq,uq, sq) be two initial data for (3.16) such that 
Uq G H m . Let U, U be the two solutions of the associated Cauchy problem. We assume 
that these solutions are defined on [0,Tq] with To > 0. Let furthermore xq G M. d , n G 
and M, Ct be as in (3.17) and (3.18). Then we have 

a (U)d t U + aj(U)djU = , 



where qq = Diag(eTo c f , 1, . . . , 1, 1) and, for all j G {1, . . . , d}, 



aj (U) = 



/ e~io<*>Uj e~<o c v ^^-(tt + \m v ) 

Uj 

■ 7D- 1 



\ 



\ 




U 3 J 



We introduce, for T G [0,Ti], 



J ^ t [(U -U)-a (U)(U -U)]+J2 d jK u -U)-aj(U)(U -U)]^ dxdt . (3.19) 



Denoting 



U-U 



(U-U)a Q (U)(U-U) and 



U-U 



U^n-Mt) (U-U)(t,x) 



dx, 
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we obtain from the Stokes formula 



I = 



[ \{U - U) ■ a (U){U - U)n t + V(£7 - U) ■ aj(U){U - U) nj I da 
JdOr \ j J 



+ 



u -u 
1 



T 







U-U 

u -u 



y/1 + 1/M 2 JC 

But we have also 

^2(U-U) aj (U)(U-U) 



M\x\ 



I loM|x| 2 v M|x| 



17-17 
1 



oM 

2 

U-U 

o 



M + e «/(2T»)c„) 2^-1(1^1 + 6|7rr)(e s/i70C ^|7r - ir\ z + \u- u\ z )^ 



, "y r */'-t">)l - 1 - I • : - ' 

(| n | +e ^/(270C„)70^1 ( | 7r | + ^| 7r |, )) 



therefore 



U-U 



U-U 



U-U 



Besides, we have 

a (U)d t (U -U) + J2 <*j{U){U -U) = J2 aoiU^aoiU^ajiU) - a Q {U)- l a 3 {U))d 3 U , 



thus 



1=1 (U -U)d t a {U)(U -u) + Y^{u -U)d jaj {U)(U -U) 
Jc T j 

+ 2 I Y.MU^ocoiUT 1 ^) - aoiU^ajWdjU . 

J Ct 



Let us denote 



h= f (U - U)d t a {U){U - U) , 

I 2 = J2(U-U)d j a j (U)(U-U), 

j 

h = 2 f Y, a o( u )(MU)~ 1 a j (U) - aoiU^ajWdjU . 
Jc T ■ 

We obtain by computing explicitly dta-o and djCtj 

h = - I —e s/ ^\TT - n)\u ■ Vs) 
Jc T 7o Ct, 
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and 



h= [ —e s/ ^ 0C "\ir-n) 2 {u-Vs)+ f div («) 

J Ct 



u-u 



7T 



(tt + &7r") V fl + (1 + 6i/7r ,, - 1 )Vir 



)) 



U-U 



dt . 



< -/i + ||DE/|| L oo(l + C(l + ||ir|| L oo)(l + ^ 

Finally, we bound I3 by computing 

aoiU^ajiU) - a {U)- l aj{U) 
I 

\ 

where VP = (e 7 o Ct> (7r + 671-") — eTo c « (jf + fevr 1 ')). We denote R = max(||s|| L<x> , 
The exponential function being convex, we have: 

< l/(^c v )e R l^\s-~s\. 



Uj - Uj 


7 V 1 (vr n + b(Tr» ^))ej T 







(Uj - Uj)I d 











Uj - Uj 



e s/(loc v ) _ gS/(70C„) 



So, 

e s/( 70 c„) (7r + _ e s/{l c v )^~ + ^ c(1 + tf-l^Wg _ 5 | + _ I) ) 

which gives us 



(l + C(l + ||7r|| L oo)(l + 6i/ 



v-l 



TT 



+ ||DC/|| LOO (1 + C(l + ||7r|| Loo )(l + 61/ 
Finally, we obtain 



7T 



)) 



»)jf 



LT — 17 



dt 



17-17 



r 



17-17 



< C(l + C(l + ||7r|| L oo)(l + 6i/ 



L°° 



)) 



dt. 



We conclude thanks to Gronwall Lemma 

_ 2 
U-U 

lB(x ,R-MT) 



(T,x)dx< \e c ' T [ 
1 Je 



B(x ,R) 



Uo-U 



(x) dx 



where C = C||DC/|| Loo(Ct) (1 + h\\^ { c T )) 
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3.2.C Estimates 

The system (3.16) can be written 

d 

A d t V + = -B(BV ,V)-J2 Cj(V)djV - 

3=1 




(3.20) 



where V = (tt, u — u, s) G M. d+2 , V = (0, u,0) and, denoting (ej)i^j^d the standard or- 
thonormal basis of R d , 

A = Diag(e s /^, 1, .... 1, (1 + t)' e ) G M d+2 
Cj = UjA , 





( 


e s /(-yoc v ) w . 


7o-i eS /(7oc„) 7r e j T 









Wjl d 







V 








(1 + t)-» Wj 



B 



( 2!i r ie s /^» c -)vrdivn \ 
[w ■ V)u 




We also introduce N k (t) = ^ J Rd B k V • A {V) B k V dx^^ and Z{t) = J2T=o( l + t ) 9kN k{t) 
with gk = k + r — a and r = 9/2 — d/2, 9 G]0, min(l, In order to obtain energy 

estimates, we apply D fc to (3.20) and we multiply it by D fc V. Then, we integrate on M. d . 
The additional term with respect to the Perfect Polytropic gases considered by M. Grassin 
[7] is now 

F* (BV , BV , V) = (div (w + u), Vvr, 0) G R d+2 . 

With the notations U = (n,w) G R d+l and F as in (3.9), the last component of F* being 
0, we have 



/ 



d K Vd k (F*(BV ,BV ,V))dx 



d k Ud k {F{BU ,BU ,U))dx 



Using the estimates (3.12)-(3.13)-(3.14) we finally get an estimate on Y^. The defini- 
tion of the norm is slightly changed with respect to isentropical case, however, for all v = 
(v 1 ...,v d+1 ) G R d+1 ,ifv* = (vi,...,Vd+i,z) G R d+2 , we have ||w|| 2 < e P"^ / A v* . 

Consequently Yj. ^ e^ s ° Hl°°/(7oc„) jy fe and the estimates on Yj, obtained in the isentropical 
case give an estimate on in the general case. 

Finally, we obtain, adding the estimate on F obtained in the isentropical case to the 
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estimates from M. Grassin in the general case: 

1 dNl + k±Z N 2 ^ c(1 + t) P N 2 z + c / (1 + t) -9u-2 NkZ + C(1 + t) /3+(/3 + l)(,-l )Ar 2^ 

— ■ cli 1 — I - £ 

+ C(l + t)(/3+i)(^-i)-iiV fe 2 ^- 1 + C(l + t)^+ 1 )("- 1 )-fffc- 2 7V fc Z ,/ 
+ C J] iV^ 2+f (1 + t )-^+/3+5(/3+i+e/2) 

where /3 = — g\ — i and 



Then, we choose a so that j3 = 0, i.e. 



a = 1 + 6/2 > 1 with e]0, min(l, -^—)] ■ 



Next, we simplify by Nk, we multiply by (1 + t) 9k and we summate on k to obtain 

^ + < c (z 2 + j^-^i + (i + ty- l z» +l + z 2+m (i + t) am ^j . 

We denote now ((t) = (l + t) a exp (iq^) Z(t) and we deduce from the inequality just above 

— ° (C 2 + C +1 + C l+2 )- 

dt " (l + t)° ^ ^ ; 

We conclude in the same way we did in the isentropical case, replacing v by v* = max(zv, m+ 
1) > 2, since ( 2 + C +1 + ( m+2 2(( 2 + 

4 Technical tools 
4.1 Lemmas 3.4 and 3.5 

We show here Lemma 3.4, which states, with the notations introduced in Section 3.1: 



|/| < CNI^HDE/Hloo 



B k U 



2 

L 2 



Proof of Lemma 3-4- 

If k = 0, I = —b f Rd tt u (-k div to + to • Vio) dx. By integration by parts, we obtain 
I = —b / it u+l div w dx ^ C tt u ~ 1 1 1 7r 1 1 ? 2 1 1 div it; 1 1 T 00 • 

If k ^ 1, I is such that 

I = -b Yl I d a 7rd Q (7^ div w) + ^ & x w j &*(it v d j ir) . 
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Expanding, we find 

I = -b Yl f d a (-ir)TT u d a (divw) + ^2d a (w j )iT»d a {d j Tr) + Z, 

aemWa^k Rd 3 

where £ is a sum of terms as J^ d d a °Ud ai (Tr u )d a2 U , where U is any of its component 
and \ao\ 1 = k, \a±\ 1 = I, \a2\1 = k+l — l, with I € [1, kj so that the derivatives are of 
order less than k. We treat first one of the terms of the preceding sum for a d-uplet 
a £ N d , of size k, that is to say \a\ 1 = a\ + . . . + = k. By integration by parts, we 
find 

d a (ir)ir u d a (dww) + ^d a {w j )it u d a {d j it) 

3 

i?Y i d j {ff*Wff*(w j )) 
3 

A J*- d 



Its. 



Hence, / is a sum of terms as / d k Ud l {Ti u )d k ' l+1 U where 1 < Z < fc - 1. 

Note that we used here the notation after which d k U means d a U for a given a G N d 
such that \a\ 1 = Yl a j = Furthermore, we write £/ for any of its component. 
Consequently, d k U can mean d a Wj. We will use by now this notation. 

If k ^ 1 and I ^ 1, I k, using the Holder inequality and the Gagliardo-Nirenberg 
Lemma 4.1, we have: 



/ 



d k ud l (ir u )d k ~ l+1 udx ^ 



d k U 
d k U 
d k U 



L 2 



L 2 



10(0 



„fc-l 




C7 


L T=T- 






Il°° 








L 2 









X IWIloo*" 1 



<9 fe ?7 



fc-i 

L 2 



We use next Lemma 4.7 and the inequality 



ll^)|| 



V1T 



dir 



< cikllL-ll^l 



L°° ' 



to obtain 



d'W(7r ,/ )d fc - z+1 i7dx< C 



L 2 



listen 



11 n^^ 1 



If k > 1 and I = k, we have to estimate 



/ 



d k Ud k (ir u )dU ^ 



d k U 
d k U 



L 2 



<9 fe K) Jdu\\ LO 



L 2 



7T 



7T 



L 2 



IWIl 
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If k ^ 1 and / = 1, we have to estimate 

f d k ud{Ti u )d k u < 



The inequality ||3 1 (tt , ')|| Loo < C , ||7r||L» 1 ||a 1 7r| 



d k U 








L 2 




2 




d k U 






L 2 



\dU\ 



L°° 



allows us to conclude. 



□ 



We prove now Lemma 3.5, which states, with the notations introduced in Section 3.1 



J 2 < C(l + t) d " 



B k U 



L 2 



where d k = {(3 + l){v - 1) - g k - 2 and /? = -51 - f . 

Proof of Lemma 3.5. For A; = 0, J2 = 0. We are considering here k ^ 1 ; J2 is then 
a sum of terms J d k Ud l (Tr")d k ~ l+1 u dx for ^ I ^ A; — 1. The choice of initial conditions 
gives us U € H m , but we do not know if D 2 U € L°°. So we can distinguish two cases: 
m > 2 + d/2 and D 2 U € L°°, or m ^ 2 + d/2. 

Case m > 2 + d/2. We study now different cases, after the values of k and I. 
1. If /c ^ 1 and / = 0, we use Proposition 3.1 and Lemma 4.7 to obtain the estimate 



/ 



TT V d k TTd k+l Tidx < |k"| 



L°° 



d k U 



L 2 



< C(l + t)^ +1 ^^ 

< C(l + t) dk 



d k+1 u 



D k U 



L 2 



L 2 



(1 + t )«*/2-fc 



-k-2 



D k U 



L 2 



2", 



since k — d/2 + 2 = /3+1+ g k + 2 we are done. 
2. If A; ^ 2 and / = 1, we have to estimate 

n u - 1 dird k wd k udx £ 



J 

JR<- 



d k ir 








L 2 





ttIIloo 1 !!^!!^ 



L 2 



< C(l + i) d/2_fe_1 (i + t f+^-^z" 



d k U 



L 2 



which concludes that case, since d/2 — k — 1+/3+ (/3+l)(z/ — 1) = — 1)— g k — 2. 

3. If k = 3 and / = 2, 



/ a 3 7r9 2 (7r^)a 2 ndx < 
JK. d 



7T 



L 2 



a 2 n 



L 2 



^ (1 + t)-3+G9+l)(f-l)-ffa 



L 2 



and we are done as 52 + 3 = 53 + 2. 
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4. If fe > 3 and I € [2, fe - 1], we use Lemma 4.2. Denoting g = 2^| and </ = 2 fc fc ; 3 1 
so that 1/g + l/(/ = 1/2, we obtain: 



/ d k 7rd l (7r u )d k ~ l+1 udx 
d l (ir u 



Since 



d k U 






L 2 


d k U 






L 2 


D 2 u 






L°° 



D 2 (tt") 



1-2/9 



D^ 1 (tt") 



2/9 
L 2 



D 2 TZ 



1-2/9' 



2/9' 
L 2 



^C(l + t)~ 3 , 



L 2 



< C{1 + t) d ' 2 - k , 



and 



D 2 7T 



L°° 



we obtain 



where 



< C(l + t)W +1 >- 2 Z" 

ill, L2 

<; C'(l + t) (,3+1)( ^ 1) " 9fc+1 Z iy . 

j 2 < c{\ + t) m A 



D ft 7T 



L 2 



m k = -3[l--) --(l+/3 + g k ) + -(((3 + l)v-2) 
qj q q 



+ (03 + l)(i'-l)-<fc + l) 

= (/3 + l)(i/-l)-^ fc -2 
= 4 • 



Case m ^ 2 + d/2. 

First, we note that the computations of the case k ^ 1 and I = 0, k ^ 2 and 2 = 1, fe = 3 
and / = 2 are similar. There remains to treat the case k > 3 and 2 ^ I ^ fe — 1. Since 
k ^ m ^ 2 + d/2 we have necessarily fe ^ 3, if d = 1, 2 or 3, and we are done. 
We assume now d ^ 4. Let us denote h = ±(fe + 1 + d/2) > 2. Then we have h ^ m and 

, , d- 1 
0<h-l^ 



1 , „ „ d-1 

2 </»-(fc + l-0<— ■ 



2 ' 2 

We introduce /&i = /i - I, h 2 = h - (fe - / + 1) and 1/qi = 1/2 - hi/d, l/q 2 = 1/2 - h 2 /d. 
Therefore 1/qi + l/q 2 = 2, which allows to use Holder inequality 

j 2 (k,i)= f d k Tid l {Tt u )d k - l+1 u&x < 



fc 7T 














L 2 






L92 
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Next, we apply Lemma 4.3 to find 



J 2 (k,l) ^ C 



8 k 7T 



L 2 



|D n (ir L 



lL 2 



L 2 



Finally, we use Lemma 4.7 and Proposition 3.1 



J 2 (M) < c 



(1 + t }d/2-h-l-g h +(P+l)(v-l) _ 



As d/2 - /i - 1 - + (/? + - 1) = (P + l){v -l)-g k -2, we are done. 



□ 



4.2 The Gagliardo-Nirenberg inequality and its consequences 
4. 2. A The Gagliardo-Nirenberg inequality 

Lemma 4.1 (Gagliardo-Nirenberg). (See [23, Prop. 3.5, p. 4[) Let r > 0, i G [0, r] anc? 
z G (L°° n H r )(M d ). Then d { z G L 2r /*(M d ) and 



d l z 



< r iirii 1 ~ i/r nr) r y\ 

^t,r || * || j_,oo 1 1 J — ^Hl2 • 



ii/r 



L 2r/i 



we /icwe D l z G L 9 /or g = 2^Hr and 



We deduce easily from Lemma 4.1 the following result. 

Lemma 4.2. Lei z € H m 6e swc/t D 2 zG L°° ; then for all k G [4, m], for all i G [2, fe], 

i fc-3 
'i-2 

< C D 2 z D* -1 * 

L9 L°° L 2 

Thanks to the Sobolev imbedding (see [23, p. 4]) and Lemma 4.1 we also prove 

Lemma 4.3. Let I g]0, d/2[ and 1/q = 1/2 — t/d. There exists C > depending on £,q,d 
such that for all z G H e (M. d ) we have 



L 2 



Proof. The space H £ (M d ) is endowed with the norm ||-|| L 2 + 
imbedding between and L 9 can the be written, for a given C > 0, 



L 2 



klli* < c(I^IIl 2 + 



D l z 



L 2 



for all z G . 



The Sobolev 



(4.1) 



Let us define now, for z G H , A G R+, the function z\ G lT such that z\(x) = z(Xx). 
Applying (4.1) to z\, and noting that 



I z a|Il<? 



A i \\z\ 



Li ' 



&z x 



L 2 



we obtain 



|z|| l ^ca^ 2 (INIl 2 + ^ 



B £ z 



L 2 



r/z 



L 2 



(4.2) 



where, by definition | — ^ = — Consequently, introducing A 



« T.2 



L 2 



, we have 



B l z 



L 2 



□ 
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Similarly, 

Lemma 4.4. Let p > d/2. There exists C > such that for all z G H p 

l 



< C\\z\ 



L 2 



2p I 



I 2p 

II 2 



Proof. We use now the continuous imbedding W 
C > such that 



C L c 



Thus, there exists 



< C(\\z\\ L2 + ||D p z\\ L2 ) , for all z G H p . 

Applying the inequality to z>, : x >-» z(Ax), we obtain, since satisfies 
and (4.2), 



N| Loo ^CA^(||z|| L2 +A p ||D p z|| L2 ). 



(II II \ i/p 
|p p zj] 2 2 ) ' we nave finished the proof. 



□ 



4.2.B Estimates 

Lemma 4.5. Lei m > l + d/2, U £ H m (R d ), r,a£R andZ be the norm defined by (3.10) 



z(t) = J2( 1 + t ) 9k v k u(t) 



k=0 



L 2 



with gt = k + r — a. Then we have: 

1. \\U(t)\\ Lac ^C(l + t)^Z(t), 

2. ||DC/(t)|| LOO ^C{\ + tfZ{t), 

3. If m> 2 + d/2, then D 2 U(t) ^ C{1 + tf~ l Z{t), 

L°° 

with (3 = —gi — d/2. 

Proof. 1. Applying Lemma 4.4 to U, we obtain: 

II^IIl- ^cII^IIl 2 ^!! " 1 ^!!^- 

But we also have 

II^IIl 2 < {i + tr 9 °z, 

so 



\B m U\\ L2 < (l + t)-*»Z, 



||^|| L oo < c(i + tyz, 

where s = - 5o (l ~ ^l) ~ (do + m)^ = -g - ^ = /3 + 1. 
2. Applying Lemma 4.4 to DU, with p = m — 1, we have 



l- 



|D?7|| L oc < C||DL r || L2 2(m - 1) IID™^!!^^ 17 
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In the same way as before, we obtain 

||D£/|| LOO <C(l + t)*Z 

where 

d d 
s = -si (1 - tt< tt) ~ (fi + m ~ l); 



2(m-l) / vyi y 2(m-l) 
= -0i"2=/3- 

3. Applying Lemma 4.4 to D 2 f/ with p = m — 2, which is possible since m — 2 > d/2, 
we finally prove the third inequality. □ 



Lemma 4.6. Lei f,<p& H m n L°°(R rf ) ; Zet a G N suc/i i/tai |a| = jfe < m. T/ien 



|^(/^)|| L2 ^C(||/|| 



L 2 



+ V 



Lemma 4.7. Lei / G (L°° n H m )(M d ). If v eN and u 2, or v e R and v e [m, +oo[, 
we have f v G H m (M d ) and, /or all a G N d suc/i i/iai |a| = ^ m, we have 



\d a n\v < cii/Hl^ 1 



D fe / 



L 2 



where C > is a constant independent from f , a, v. 



Proof. For f G N, z/ ^ 2, we proceed by iteration on v, using Lemma 4.6. 
For v G [m, +oo[, we have 

d a (r(x))= J2 E c a ,pr-id^f...ahf. 

l^j^k /3i + ...+/3.,-=c« 

a 2 — 

Then we take d Pi f G L & i , we apply Holder inequality and Lemma 4.1 to obtain: 



(4.3) 



3=1 



v-3 



j=i 



n k/ 

i=l 
j 

n 



'57 



^I/IIl- 



D fe / 



i=l 



L 2 



L 2 



using besides that for all j G [l,m]], 



< II/IIl- since ^ - j ^ o. 



□ 
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